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Critical collapse in tensor-multi-scalar gravity theories is studied, and found that for any given target 
space all the theories conformally related belong to the same universal class. When only one scalar 
field is present, the universality is extended to include a class of non- linear gravity theories. 
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Studies of black hole formation from gravitational col- 
lapse of a massless scalar field have revealed interesting 
non-linear phenomena at the threshold of black hole for- 
mation PJ. Specifically, starting with spherical space- 
times, 
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where dfl = d9 2 + sin 2 Bdip 2 , and {x^} = {t, r, 6, ip} are 
the usual spherical coordinates, Choptuik has investi- 
gated the Einstein-scalar field equations and found the 
following striking features: Consider a generic smooth 
one-parameter family of initial data, say, S[p], such that 
for p > p* , black holes are formed, and for p < p* no 
black holes are formed. Then, (a) the critical solution 
is universal with respect to the families of initial data 
considered, and periodic, A*(^, t) — A*(£, t + A), where 
A* = {a*,a*,(f>*}, and 
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with ro being a dimensionful constant, <f> the massless 
scalar field, and £o( T ) a periodic otherwise arbitrary func- 
tion with period A. The constant A is a dimension- 
less universal constant, which has been numerically de- 
termined as A w 3.447. (b) Near the critical solution 
but with p > p* , the black hole masses take a scaling 
form Mbh = ^(p — p*) 7 , where K is a family-dependent 
constant, but 7 is another dimensionless universal con- 
stant with 7 w 0.37. The above results were soon con- 
firmed by several independent studies both numerically 
and analytically ||, and extended to other matter 
fields, such as, axisymmetric gravitational waves per- 
fect fluids ||, non-linear er-models HQ, Einstein- Yang- 
Millers fields ||, Einstein-Maxwell-scalar fields |J, and 
more recently Einstein-Skyrme fields JlO| . 

Now it is clear that in general the critical solution and 
the two dimensionless constants A and 7 are universal 
only with respect to the same matter field, and usually 
are matter-dependent. Moreover, the critical solutions 
can have discrete self-similarity (DSS) or continuous self- 



similarity (CSS), or none of them, depending on the mat- 
ter fields and regions of the initial data spaces JTl| . So 
far, in all the cases where the critical solution either has 
DSS or CSS, the formation of black holes always turns on 
with zero mass, the so-called Type II collapse, while in 
the cases in which the critical solution has neither DSS 
nor CSS, the formation always turns on with a mass gap, 
the so-called Type I collapse. We shall first consider type 
II collapse, and then turn to type I collapse. The uni- 
versality of the critical solution and the exponent 7 now 
are well understood in terms of perturbations JL2| , while 
the one of A still remains somewhat of a mystery. The 
former is closely related to the fact that the critical so- 
lution has only one unstable mode. This property now 
is considered as the main criterion for a solution to be 
critical Q. 

In this Letter, we shall consider the critical collapse 
in tensor-multi-scalar (TMS) and non-linear (NI) grav- 
ity theories, and pay special attention on universality of 
the collapse, that is, systems that have the same critical 
behavior near their critical points. This is particularly 
interesting in gravitational collapse, as the field equa- 
tions involved are usually very complicated Jl3| , and if we 
could find a universal class, then we only need to study 
the simple ones of the class in order to know the proper- 
ties of others that usually are too complicated to be stud- 
ied either analytically or numerically. Yet, the TMS and 
NI gravity theories are the most promising alternatives 
to Einstein's theory, and are well motivated from super- 
strings and quantum gravity 14 1. It has been shown that 
the TMS theories can equally pass all the known obser- 
vational and experimental constraints as GR does, while 
can be significantly different from GR in the strong field 
regime. These differences are within the scope of the un- 
dertaking experiments, such as, LIGA and VIRGA |jl5| . 
NI gravity theories are usually connected with quantum 
corrections and are essentially also related to the strong 
field regime fig] . Since critical phenomena are strong 
field ones, it would be very interesting to study critical 
collapse both in TMS and NI theories. 

So far, only two universal classes have been found, one 
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is the massless scalar field with any potential |jj , and the 
other is a fluid with the equation of state, p = kp + f(p) 
for any smooth function f(p) Jl2| . However, they are 
all within the framework of GR. In this Letter we shall 
provide the third class but in terms of different theories 
of gravity. 

The held equations for the TMS theories are given by 
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where cj) a (a = 1,2, ...,N) denote the N scalar helds, 
B(4>) the potential, B a (4> c ) = dB{4> c )/d^) a , and j ab (j>) 
is the metric in the target space, while J^ c (4>) are the 
corresponding Christoffcl symbols. From now on we shall 
use quantities with bars to denote the ones defined in the 
Einstein frame, and quantities without bars to denote the 
ones defined in the Jordan frame which will be considered 
as the physical one In the following, we shall 

first consider the cases N = 1, 2, and then generalize the 
obtaining results to the case for any N. 

When only one scalar field is present, Eqs.(|J) and 
(Q) reduce to the corresponding Einstein-scalar equations 
considered in making the conformal transformation 
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we find that S\g,<f>] = J(~g) 1/2 (R - <j> , Q <K)dx 4 = 
J(-g) 1/2 [f(<t>)R ~ h{(j))(j). a (jy a ]dx i , which shows that in 
the Jordan frame the action represents the general 
Bergman- Wagoner-Nordtvedt (BWN) tensor-scalar the- 
ories flilJilill , where /(</>) and h((f>) are strictly positive and 
non-singular, and a prime denotes the ordinary differen- 
tiation. The BWN theories include several interesting 
cases, such as, /(</>) = 1 — A0 2 , h{4>) — 1. This is the 
only case where the coupling constant A is dimensionless. 
In this case critical collapse has been studied numerically 
jl7j and found that the two constants A and 7 are weakly 
dependent on A, while by dimensional analysis some kind 
of dependence has been expected Another inter- 

esting case is f{4>) — <f), h((f>) — u>/<f>, which corresponds 
to the Brans-Dicke (BD) theory, where u> is the BD cou- 
pling constant. In this case critical collapse has been also 
studied numerically jis|] , and found that A depends on ui 
weakly while 7 strongly. We shall show that all the BWN 
theories belong to the same universal class and have the 
same A and 7. 

To show the claim, let us first prove that the solution 
A* = {<?„„, 0*}, which is obtained through Eq.(||) from 
the critical solution A* found in 0], is a critical solution 
in the BWN theories. First, since f(4>) is non-singular, 



the solution A* has the same singular behavior as A* 
according to Penrose's theorems |19|| . In particular, if 
A* is no singular on the axis r = for any t < 0, so is 
A*. Moreover, ^ is a monotonically increasing function 
of <j) [otherwise, we can restrict the mapping to a specific 
region, so that it is one-to-one], so the inverse (j> = <t>{4>) 
exists and <f>* has the same period as </>*. Thus, if A* has 
period A so does A*. Second, if the perturbations of the 
critical solution A* are given by || 
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where 5iZ(^,r) are periodic functions with period A, 
then from Eq.([|) it can be shown that the perturba- 
tions of the solution A* to first order take the same 
form only with Z being replaced by Z in Eq.(|5|), where 
Z = {a,a,(f>.£,(f>. T } ||. Since the equations satisfied by 
SiZ are related with the ones satisfied by SiZ through 
the regular transformation (|J), the perturbations have 
the same spectrum in both cases. Thus, if SiZ has only 
one unstable mode, say, Ai, so does SiZ. Therefore, if A* 
is a critical solution of the Einstein-scalar field equations, 
A* is a critical solution of the BWN field equations. 

Assuming that the linear perturbations are valid until 
the moment r — To, we have 

Z p {r) = Z(£, r ) = Z*Qjx —,7ft) + eStZQn. -,r ), (6) 



where r p = roe r ° + ^ ^ r °', and e 



(p-p*)e AlT0 5Ci(p*)/5p 
is very small so that the perturbations are valid. Note 
that in writing the above expressions, we have dropped 
all the stable modes and expanded C\ (p) to its first order 
by using the fact that C\(jp*) = 0, because when p = p* 
we shall obtain exactly the critical solution. To show the 
explicit dependence on r p , we have also replaced £ by 
\n(r/r p ). Now we need to use Eq.(|^) as the new initial 
data to consider the non-linear evolution of the BWN 
field equations, which in the Jordan frame read 
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Y = 



ee (27rG) 1 / 2 ^ r /a 
(27rG) 1 / 2 0, t /a, W = f-^ 2 d^j)/d(j), X± = X ± Y, and 
C\, C2, and C± are functions of a, f, W and X± only. 
Now the key observation is that the initial data Eq. (||) de- 
pend on r p only through the form r/r p and to/r p , while 
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the field equations (]?]) are invariant under the transla- 
tion of r. Then, the entire solution must take the scale- 
invariant form flql 



Z(tT) = Z(ln-,-). 

Tp Tp 
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On the other hand, we know that the mass of black holes 
has dimension length, while in the above solution the only 
quantity that has such a dimension is r p . Therefore, we 
must have 



M BH (xr p = r e^ +T1 (p-p*) 
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where 7 = — A^ , n = e~ 1 j\n[dCi(p*) /dp}. The above 
analysis shows that the gravitational collapse of the BWN 
scalar field in the BWN theories exhibits the same crit- 
ical behavior as that of a massless scalar field in Ein- 
stein's theory. In particular, they have the same A and 
7. Therefore, they belong to the same universal class. 
Now let us turn to the NI theories given by the action 
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where f(R) is an analytic function of the Ricci scalar 
R. In this case the collapse of the gravitational field 
also exhibits the same critical behavior as the mass- 
less scalar field in Einstein's theory studied in Q|. This 
is because by the conformal transformation e x °^ — 
f'(R), g„ v = e x '^g^, we find that S non .[g,R} = 

fd*x(-g) 1/2 f(R) = I d*x(-g) 1/2 lR - fr*4>, a - 2V(4>)], 

where V((f>) = e - Ao *[i?(0) - e"^ f(R(4>))]/2, and A = 
(2/3) 1 / 2 . That is, the non-linear gravity theories given 
by the action (|l^) is related to the Einstein-scalar theory 
by a conformal transformation, where the scalar field <p 
usually has a non- vanishing potential V(<fi). As shown 
numerically in and analytically in p"2] |, the potential 
does not affect the (type II) critical behavior, so for any 
given potential the collapse will exhibit the same critical 
phenomena as the one with vanishing potential. Then, 
using the same conformal transformation arguments as 
given in the BWN case, we can show that the critical col- 
lapse of the gravitational field in the NI theories exhibits 
the same critical behavior as a scalar field with a po- 
tential in Einstein's theory. Therefore, the gravitational 
collapse of a pure gravitational field in the NI theories 
given by Eg. JTZjJ and the one of the BWN gravitational 
scalar field in the BWN theories all belong to the same 
universal class as that of a scalar field ( with or without 
potential) in Einstein's theory. 

When two scalar field are present, a particular case was 
studied in H[7|, where the functions j a b(4> c ) was chosen as 
lab = A" 1 diag.{l,g(^)}, with gffi) = e -*\ sin 2 20\ 
and A being a coupling constant. When g = , the 
non-linear a— model is equivalent to the BD theory cou- 
pled with a massless scalar field 4> 2 with A = 8k — 
(u> + 3/2) -1 > 0, which has been studied in details in 



0, while when g — sin 2 20 1 , it is the case studied in || 
for k = A/8 < 0, where n is a constant. As noticed in 
H , in the latter case there is no correspondence between 
the NI a— model and the BD theory. For more_details, 

It has been shown 



we refer readers to 20 . It has been shown |||7| that 
critical (Type II) behavior exists for all the range of k, 
although the properties of the collapse crucially depend 
on k. In particular, For k > 1/12, a family of critical so- 
lutions characterized by k exists, and all these solutions 
have CSS. The exponent 7 strongly depend on k in this 
range. For < K < 1/12, the critical solution found in 
jl| seems to be the only one, and both 7 and A are now 
weakly dependent on k. For k < —0.28 further bifurca- 
tions are present. In HQ, the problem was studied in 
the Einstein frame. An open question is that do the re- 
sults still hold in the BD theory (or in the Jordan frame) 
jjl| ? as it is well-known that the mass of black holes are 
not well defined in the BD theory [jl4|]. In the following, 
we shall use the conformal transformation arguments to 
show that this indeed is the case. 

The conformal transformation from the Einstein the- 
ory to the BD theory is given by g^ v = f(4> 1 )g^ with 
/(</r) = (ft 1 . However, to have the results applicable as 
much as possible, we shall take f{<j) 1 ) as a non-singular 
otherwise arbitrary function of </r. Moreover, we will 
not consider only the case where k > 0, but also the 
case where K < 0, although in the latter case no known 
theories are related. Since the conformal factor /(<?r) is 
non-singular, we can see that the metric in the Jordan 
frame will have the same singular behavior as that in 
the Einstein frame. Yet, if the solution Z* is a critical 
solution in the Einstein frame, the solution Z* is a crit- 
ical solution in the Jordan frame, because the equations 
for the perturbations 8Z are related to the ones for SZ 
by a regular conformal transformation, and consequently 
they have the same spectrum of the unstable modes in 
both cases, where Z represents a set of variables of fields 
such that the complete field equations can be written 
in the first-order form F(Z,Z, T ,Z^) = 0. Moreover, if 
this equation is scale-invariant under the translation of 
r, so is the equation F(Z, Z_ T , Z^) = for Z. Then, us- 
ing the scale- invariant arguments given previously we can 
see that the entire solutions must take the scale-invariant 
form Eq.(||), and as a result, the mass of black holes for 
the supercritical case will take the same scaling form in 
the both theories. 

Clearly, the above analysis does not depend on the 
particular choice of 7ob(</> 1 , 4> 2 ) and f (cj) 1 , (f> 2 ) . Therefore, 
for any given theory, that is, for any given target space 
7a&(0 > 4> 2 ), au the theories related to it by the confor- 
mal transformation g^ v — f(4> 1 ,<p 2 )g[iL> belong to the 
same universal class. Moreover, if the critical solution 
Z* has CSS (DSS), the solution Z* has CSS (DSS), too. 
In the DSS case, a "wiggle" overlaid on the scaling law 
was found || . It is not difficult to see that it will persist 
after the conformal transformation. 

It is not difficult to be convinced that the above con- 
clusions do not depend on how many scalar fields cf> a in- 
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volved and the explicit form of f((f> a ), as long as f{4> a ) 
is positive and non-singular. Also, they do not depend 
on the symmetries, such as, spherical, axial, etc., of the 
spacetimes. 

For type I collapse, using the same conformal trans- 
formation arguments as given in the type II case, it can 
be shown that if Z* is a type I critical solution in the 
Einstein frame, Z* is also a type I critical solution in the 
Jordan frame. This is because in this case the perturba- 
tions 5Z will take the form lO] 



6Z(t, r) 



dp 



e Xl %Z(t, r), 
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in the Einstein frame. Clearly after the conformal trans- 
formation, the perturbations 5Z to the first order will 
take the same form. Since the equations satisfied by SZ 
are conformally related to the ones satisfied by 5Z, they 
must have the same spectrum of the perturbation modes. 
Thus, if Z* is a critical solution in the Einstein's frame, 
Z* is a critical solution in the Jordan frame. On the other 
hand, Eq.(ll) is different from Eq.(||) in that in the for- 
mer the perturbations are proportional to e Al *, while in 
the latter they are proportional to (— i) . It is this dif- 
ference that determines the types of the collapse. Since 
the conformal transformation does not change the forms 
of perturbations, it is concluded that it does not change 
the types of collapse, and that all the theories that are 
related by a conformal transformation belong to the same 
universal class for type I collapse, too. 

In conclusion, we have shown that all the theories that 
are related by a conformal transformation belong to the 
same universal class, irrespective of the spacetime sym- 
metries. The conformal transformation preserves the 
types, Type I, or Type II, of the collapse and the symme- 
tries, CSS or DSS, of the critical solutions. As a result, 
all these theories are distinguishable near their critical 
points. 
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